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Abstract 

We discuss the possibility of the "mixed state" in incommensurate charge 

density waves with three-dimensional order. It is shown that the mixed state 

can be created by applying an electric field perpendicular to the chains. This 

state consists of topological dislocations induced by the external field and is 

therefore similar to the mixed states of superfluids (type-II superconductor or 

liquid Helium II). However, the peculiar coupling of charge density waves with 

the electric field strongly modifies the nature of the mixed state compared to 

the conventional superfluids. The field and temperature dependence of the 

properties of the mixed state are studied, and some experimental aspects are 

discussed. 
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The charge density wave (CDW) is an ordered state in which translational symmetry is 
broken [|TJ. Topological dislocations (TD's) of incommensurate CDW with three-dimensional 
(3D) order have been studied by several authors with connection to the phase slip phenomena 
in sliding CDW's [0-[J. The analogy of the TD's with the dislocations in crystals has also 
been pointed out |||7|]. However, unnoticed up to now seem to be the similarities of the 
TD's with the vortices in superffuids: both are topological singularities of complex order 
parameter. In type-II superconductors in magnetic fields and liquid Helium II in a rotating 
container, the vortices create a remarkable, so-called mixed state, which generates a wide 
variety of intriguing phenomena ||. In this Letter, we investigate the possibility of the 
corresponding mixed state in CDW's and discuss its properties. 

In contrast to the vortices in type-II superconductor and liquid Helium II, which induce 
magnetic flux and angular momentum, respectively, the TD's in CDW's induce charge po- 
larization. Therefore, in analyzing the electrostatic properties of the TD's, we have to treat 
the scalar potential carefully. We start with the microscopic model of the electron-phonon 
system in the presence of scalar potential, and derive the Ginzburg-Landau free energy of 
the ordered state. The effective free energy of the TD's is obtained by integrating out the 
single- valued part of the phase of the order parameter, i.e., the phasons. The TD's in the 
CDW's are similar to the vortices in type II superconductors or superfluids, except for the 
coupling to the external field. We show that the mixed state is created by an external 
electric field in the direction perpendicular to the chains (transverse direction), which can 
be understood as follows. When the system tries to screen the field, the chemical potential 
must be changed from chain to chain so as to induce a surface charge |J. This modulation 
changes the wave number of the CDW condensate, i.e., 2kp, thus conflicts with the 3D order 
and causes frustration. We analyze this state based on the free energy of TD's. First, we 
determined the lower critical field D c i at which the first TD appears in the system as we 
increase the field strength. If we increase the field further, the density of the TD's also 
increases. In the strong field region, the density of TD's become so high that we can treat 
them as a continuum. In this limit we find that the width of the mixed state is given by the 
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Thomas- Fermi screening length of the normal state. If we increase the field even further, 
the cores of the TD's begin to overlap, thus destroying the 3D order. The characteristic field 
of this phenomena, we denote by D c2 , is estimated. We argue that these properties should 
be experimentally observable. 

First we derive the free energy of CDW from the ID electron-phonon system. The 
imaginary-time action reads (here we act in the following units H = ks — 1 unless noted), 
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where ip a = l/vLJ2k e lkXa k,o and b q are variables expressing the electron and phonon degrees 
of freedom with excitation spectra = k 2 /2m and u q , respectively. L, g q , (3 and \x are the 
length of the chain, the electron-phonon coupling constant, the inverse of temperature, and 
chemical potential, respectively. The scalar potential is expressed by <p and the charge of an 
electron is given by — e. From the above action, the free energy of the CDW can be derived 
in both the cases of T<T C and T <ti T c . In the former case, we utilize the expansion of the 
order parameter up to fourth order and the scalar potential up to second order. On the other 
hand, in the latter case, we use a gradient expansion of the phase of the order parameters 
and the scalar potential |[LU|| . When the quasi-particle excitation gap A is uniform, the free 
energy becomes in both cases simply, 



F = — (1 - f s ) f dxip 2 + i-f s { dxOE 

+ V ^f s J dx{d x 6f (2) 

disregarding irrelevant constants. Here 9 is the phase of the order parameter, E = —d x ip 
and f s = |A| 2 7rT£e n (|A| 2 + e 2 )~ 3/2 is the condensate density given by 4c |A| 2 for T<T C 



with c = /3 2 C(3, 1/2)/(4tt) 2 , and 1 - ^2vr/5| A|e"^ A l for T < T c , where C(3, 1/2) is the 
zeta function. The first term of Eq. (0) expresses the screening due to the excitation of 



quasi-particles [fLl]] . The second term describes the acceleration of the CDW and coincides 



with the result obtained from the "chiral transformation" at T <C T c [f[2|,[n| . 

The free energy of the 3D ordered state can be obtained by introducing the rigidity due 
to the inter-chain coupling, as follows, 

K r, 
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where the electrostatic energy is also included. We assume that the chains are parallel to 
the x-axis. The anisotropy is parameterized by 7 = £ y /£ x and 7' = £ z /£ x with £ x , £ y and 
£ z being the coherence lengths in x-, y- and ^-direction, respectively. The screening length 
due to the quasi-particles, Ao, is given by Aq 2 = 8A^"j_e 2 (l — f s )/vF- The response of the 
system is probed by the external charge density p ext (r) which induces an external field given 
by D^f* = Aneikp^ 1 / k 2 where k 2 = |k| 2 . K = N_ L v F f s /2n and J = N ± f s /n are coefficients 
proportional to the areal density of the chains N±. Note that the size of the cores of the 
TD's is given by £ y and (, z in x-, y- and z-direction, respectively. Although the spatial 
variation of the amplitude of the order parameter is neglected, the present treatment is 
applicable to the most CDW systems for the following reason. Since the transverse size of 
the cores is usually smaller than the inter-chain spacing except near T c , the dislocations 
mostly sit between the chains so as to minimize the free energy, and the cores thus do not 
affect the order of CDW. In this Letter we focus on the intrinsic properties of a clean CDW 
condensate and leave the effects of inhomogeneity pinning for further work. 

Based on this free energy, we first clarify the screening properties of the CDW state 
without TD's. In this case, only the phason modes contribute to the screening and leads 
the free energy: 

F eff =Y\jrz (k 2 + A k 2 ) Irf + iew^l , (4) 

where the effective screening length A k is given as Aj^ 2 = Aq 2 + (Ane 2 J 2 /K)k 2 /k 2 and 
k 2 = k 2 + l 2 k 2 + 7 ' 2 /c 2 . In the low temperature limit, Ao diverges since the quasi- 
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particle excitations are exponentially suppressed by the energy gap, and then only the 
polarization of the condensate can contribute to the screening. In the direction along the 
chains (longitudinal direction: k y = k z = 0), the phason contribution completely compen- 
sates the suppressed quasi-particle contribution. Actually, the screening length becomes 
(Aq 2 + 47re 2 J 2 /i^) -1 / 2 , which coincides with the Thomas-Fermi length of the normal state, 
X TF = {8N±e 2 / (hvp)}^ 1 ^ 2 . On the other hand, in the transverse direction [k x = 0), there is 
no phason contribution and, consequently, no screening for T <C T c . This can be attributed 
to the "rigidity" of CDW. 

Next we examine how the TD's affect the screening properties. Here we limit our 
discussion to straight dislocations parallel to the z-axis, which reduces the problem to a 
two-dimensional one. The density of topological charge is then n(r) = Y.fj,Q^ 2 K r ~ r n), 
where g M and are the topological charge, = 2n x (integer), and the position of the \i- 
th dislocation, respectively. The phase of the order parameter is given by the relation, 
(V0)k = «k#£ + ik x z n^/k 2 , where 0£ expresses the phason part, z = (0,0,1), and 
k = (k x ,k y ,0). In the following we consider TD's with = ±2ir only. By integrating 
out the phason mode, the free energy of TD's per unit length becomes, 




where L z is the length of the sample in the ^-direction. Note that in the imaginary time 
formulation the scalar potential has to be rotated to the imaginary axis like y?k — > icp^. 
Therefore we define the expectation value of ip^ by (y?k) = ^k, where <^k is the physical 
value. 

Eq. (|5|) tells us that the TD's behave like a Coulomb gas with background charge 
—eJ(k y ip\t) I K , which is proportional to the electric field in ^/-direction. Therefore TD's are 
generated by the external electric field in ^/-direction, forming the "mixed state" of CDW. 
In order to quantify this point, we study two limiting cases: the dilute and the dense limit 
of dislocations. 
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In the dilute limit, we consider the problem of a single TD: When does the first TD 
appear as we increase the field strength? The corresponding critical field can be estimated 
from the difference of the free energy between the single TD state and the TD free state, 
AF. It is given in the limits of T <C T c and T<T C as, 



AF kK^WXtf - 27ceJD ext W 2 for T < T c , 

(6) 

nK-y ln(W/fv) - 2neJX 2 TF Df for T<T C . 
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In calculating Eq.(||), we assumed that the system is infinite in the x-direction and of width 
W in the ^-direction, thus introducing an infrared cut-off of 1/W for the wave vector k y . Note 
that "T<T C " should not include the region too close to T c , where the present treatment fails 
due to the diverging coherence length. In the limit of T <C T c , the critical field is expressed as 



D c \ = ( / -f/4e)hu) p xW~ 1 , where uj p is the plasma frequency given by \J8N±e 2 VF/h. Therefore, 
the critical voltage defined by V c \ = D C \W becomes independent of the sample width. On 
the other hand, at T<T C , D cl is given by / Ae)hoj p x X^p \n(W/ Xtf)- This is larger than 
the D c \ of the low temperature limit, because the energy gain due to the creation of TD's 
is reduced by quasi-particles. 

In the dense limit, n-^ is approximated by a continuous function rather than a function of 
the positions of individual TD's {r M }. This corresponds to neglecting the spatial structure 
smaller than the spacing between TD's. In this approximation, it is quite simple to evaluate 
nk and (p-^ by using Eq. (||). The minimization of the free energy F e ff in terms of gives, 

n k = -^kytpk. (7) 

The density of TD's turns out to be proportional to the strength of the electric field in the 
y-direction. On the other hand, the variation of F e ff with respect to (py_ gives the Poisson 
equation, 



(k 2 + A k 2 )^ k = -47re ( J^ikyU^ + ) . (8) 
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Substituting Eq. (0) into Eq. (g), we obtain, ip k = —Anep^/ik 2 + A m 2 ), where 
Aq 2 + Aire 2 J 2 / K . Note that A m equals now the Thomas- Fermi screening length of the 



normal state, Xtf- When we apply the electric field in the ^-direction, the field penetrates 
the sample to the depth A m . As we can see from Eq. ([?]), the TD's are thereby distributed 
in this surface region where the mixed state is formed. 

In case of type-II superconductors (or liquid Helium II), there is another characteristic 
field, i.e., the upper critical field (or angular velocity), usually denoted as H c2 (fi c2 ) [§, 
which defines a phase transition within the mean field theory. H C 2 (fi C 2) can be understood 
as the field at which the cores of the vortices begin to overlap. In case of the CDW's, 
when the cores of the TD's begin to overlap, the inter-chain ordering is almost destroyed 
by strong fluctuations due to the dimensional reduction. We denote this characteristic field 
as D c2 for convenience, although it is still unclear whether there is a real phase transition 
or not. This field strength is first realized at the surface where the internal field coincides 
with the external one. Since the cores begin to overlap when the number density of the 
TD's exceeds l/(£ x £,y), the corresponding field strength, D c2 , is estimated from Eq. (0) as 
D C 2 = nv F h/(e£ x Q. 

In the following, we discuss the experimental aspects of the mixed state with actual 
CDW materials in mind. We concentrate on the low temperature limit where the effects 
due to the mixed state are rather clear. Since the TD's modify the screening properties of 
the system, the mixed state can be observed by measuring the dielectric constant of the 
system in a capacitor with the transverse electric field. The capacitance should display the 
crossover from the insulating behavior in the weak field region (V < V c i) to the metallic one 
in the strong field region (V > V c i). In the case of K 3 MoO 3 (7 = 0.01 ~ 0.1 depending 
on the direction U) , V c \ is estimated as 7 x 0.69 V. On the other hand, D c2 is estimated 
as 7 _1 x 2.9 x 10 6 V/m, which can be much larger than V& even if we take W = 1 fim. It 
should be noted that, in the case of K 03 MoO 3 , a simple estimation of the Thomas- Fermi 
screening length in the normal state, Xtf — ^p/vf, gives less than 1 A, which is smaller 
than the inter-chain spacing ~ 10 A. This estimation is based on the continuum model of 
metals and, therefore, we cannot take it too seriously since the continuum approximation 
is no longer valid at this scale, especially in the transverse direction. However, we expect 
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that, in the strong field region, the width of the mixed state can be as narrow as the lattice 
spacing and most TD's are confined to the outermost layers. In the weaker field region, this 
width can be much wider. 

In the uniform electric field, there can be an additional effect due to the quasi-particles 
excited through Zener tunneling process |TJ]]. In the WKB approximation [15], the char- 
acteristic electric field of Zener breaking is huge in the present configuration (~ 7 _1 x 10 7 
V/m) due to the narrow band width and the large lattice spacing. Therefore the contri- 
butions from the quasi-particle excitations is strongly suppressed and we expect that the 
mixed state of the condensate can be observed. 



The effects of pinning cannot be avoided in realistic CDW systems |T(|, but are beyond 
the scope of the present Letter. However, we expect a hysteretic behavior in the observed ca- 
pacitance. Of course, the magnitude of the hysteresis depends strongly on the inhomogeneity 
of the sample. 

In conclusion, we proposed a new state, i.e., the mixed state, in CDW's, which is re- 
alized by applying the transverse electric field. The physical properties are quantitatively 
discussed based on the Ginzburg-Landau free energy of the CDW and the critical fields which 
characterize the mixed state are estimated. The mixed state strongly affects the screening 
properties of the system and is thus detectable through the dielectric constant measurement 
in a capacitor. 

The authors wishes to thank Prof. H. Fukuyama and Prof. H. Matsukawa for stimulating 
discussions and, Prof. G. Bauer, Dr. Y. Nazarov and Dr. C. Dekker for useful comments. 
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